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ABSTRACT 

Wc measure the acoustic scale from the angular power spectra of the Sloan Digital Sky Survey HI 
(SDSS-III) Data Release 8 imaging catalog that includes 872, 921 galaxies over ^ 10, OOOdeg^ between 
0.45 < z < 0.65. The extensive spectroscopic training set of the Baryon Oscillation Spectroscopic 
Survey (BOSS) luminous galaxies allows precise estimates of the true redshift distributions of galaxies 
in our imaging catalog. Utilizing the redshift distribution information, we build templates and fit to 
the power spectra of the data, which are measured in our companion paper, iHo et al.l (|2011[ ). to derive 
the location of Baryon acoustic oscillations (BAO) while marginalizing over many free parameters to 
exclude nearly all of the non-BAO signal. We derive the ratio of the angular diameter distance to the 
sound horizon scale D^(z)/rs = 9.212to'4|!)4 at z = 0.54, and therefore, Da{z) — 1411±65 Mpc at z = 
0.54; the result is fairly independent of assumptions on the underlying cosmology. Our measurement 
of angular diameter d istance Da{z) is lAa higher than what is expected for the concordance ACDM 
(|Komatsu et al.ir2Qllt) . in accordance to the trend of other spectroscopic BAO measurements for z > 
0.35. We report constraints on cosmological parameters from our measurement in combination with 
the WMAP7 data and the previo us spectroscopic BAO measurement s of SDSS (iPercival et al.ll2010D 
and WiggleZ (|Blake et al.ll2011b[ ). We refer to our companion papers (|Ho et al.ll2011l : lde Putter et al.l 
120 lit ) for investigations on information of the full power spectrum. 
Subject headings: cosmology 
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1. INTRODUCTION 

Baryon acoustic oscillations (BAO) imprint a distinct 
feature in the clustering of photons (i.e., cosmic mi- 
crowave background), mass, and galaxies. Sound waves 
that propagated through the hot plasma of photons and 
baryons in early Universe freeze out as photons and 
baryons decouple and leave a characteristic oscillatory 
feature in Fourier space and a single distinct peak in 
the correlation function approximately^^ at the distance 
the sound waves have traveled before the epoch of re- 
combination. The distance is called the "sound horizon 
scale" and determines the physical l ocation of the BAO 
feature in clustering statisti c s (e.g.. [Peebles fc Yul 11970 ; 



Sunvaev fc ZeldovichI [19701: IBond fc Efstathiou 



1984: 



HoltzmanI 119891: IHu fc Sugivamal Il996t IHu fc Whitg 



MoltzmanI I19»a IMu & Sup 
1996 HEisenstein fc Hulll998ir 
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The scale observed in the mass is not exactly the distance trav- 
elled when recombination occurs as the momentum of the baryonic 
material means that the motion continues for a short time after 
recombination, until an epoch known as the baryon-drag epoch. 
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Cosmic microwave background (CMB) data provides 
an independent and precise determination of the sound 
horizon scale. Therefore, comparing this sound horizon 
scale to the observed location of the BAO from galaxy 
clustering statistics allows one to constrain the angular 
diameter distance and Hubble parameters, thereby 
providing information on the nature of dark energy. 
This approac h is known as the 'standard ruler test' (e.g., 
Hu fc Whit^ p96: Eiscnstcin 2003; Blak e fc Glazebrook 
2003|lLindeJ[2003tlHu fc Haimanii2003t iSeo fc Eisenstein 
2003| ). BAO techniqu e is considered an esp ecially robust 
dark energy probe (jAlbrecht et al.l I2006f ) for various 
reasons. First, its physical scale is separately measured 
from CMB data. Second, the nonlinear effects in the 
matter density field are still mild at the BAO scale 
(~ 150 Mpc) such that the resulting systematic effects 
are small and can be mo deled with low-order per- 
turba tio n theories (e.g. , iMeiksin. White, fc Peaco"ck 
1999t ISeo fc EisensteinI 120051: iJeong fc Komatsu 
2006 : iCrocce &: Scoccimarrol l2006t lEisenstein et al. 
2007 : iNishimichi et al.l l2007t ICrocce fc Scoccimarrc 



20081 iMatsubara 
Seo et all l2nll8P 



200; 



7. ' 



Padmanabhan fc White! I200£ 
Taruva ct al. 2009: Seo et al.l 12011" 



Third, the observational/astrophysical effects such 
as galaxy /halo bias and redshift distortions are 
likely smooth in wavenumber and do not mimic 
BAO s uch that they can be marginalized over (c.g^ 
Seo fc Eiscnstci n 2005; Huf f et al. 2007; Sanchez etaL| 
20081 IPadman abhan fc Whitd 12009; Mehta ct al. aom 
[but see Dalai et al.l ()2010f ) and lYoo et al.. (2011. ) for a 
possibility of an exotic galaxy bias effect]. 

In recent years, BAO have been detected in the 
galaxy d istribution and used t o constrai n cos- 



Hatsil 120061: iTeemark et al.l 120061: iPercival et al. 


2007allbl; Padmanabhan et al.l 20071; iBlake et al. 


20071 lOkumura et al.l 120081; lEstrada et al.l 


2009: 


Gaztanaea et al.l l2009allbl: IPercival et al.l 


2010: 


Kazin et al. 201C 


: Beutler et al.l 120111: ICrocce et al. 


201lHBlake et al.l 


2011al|br). Most of these studies have 



used a 3D distribution of galaxies from spectroscopic 
surveys to constrain an isotropic distance scale Dv{z) 
{Dv{z) = [(1 + zfD\{z)cz/H{z)Yl^ where Da is 
the angular diameter distance and H is the Hubble 
parameter) using spherically averaged clustering statis- 
tics, while others have constrained Da{z) and H{z) 
separately, using anisotropic clustering information. 

Retrieving 3D spatial information requires accurate 
redshift determination (i.e., spectroscopic surveys), de- 
manding specialized spectrographs and surveys that typ- 
ically take longer times. Multiband imaging surveys, 
on the other hand, can more quickly cover a large 
number of galaxies (low shot noise) and a large area 
of sky but provide only 2D spatial information, as- 
suming a realistic level of photometric redshift error, 
and therefore fail to retrieve information on iJ(z)^^. 
Another disadvantage of using the imaging data to 
make BAO measurements is an additional damping of 
the BAO due to projection effects and difficulty in 
applying BAO reconstruction. Nevertheless, imaging 
surveys can, in principle, provide larger and deeper 

A fractional redshift error of 0.25% in l-fz is at least required 
to recover H{z) llSeo fc Eisensteinll2003l) . 



surveys (ISeo fc EisensteinI l2003t lAmendola et al.l 120051: 
Blake fc Bridld 120051; iDolngy et al. I I2006I; IZhan fc Knoxl 
20060 and this prospect has motivated current and fu- 



ture imaging B AO surveys such as the Dark Energy 
Survey^^ (DES; iThe Dark Energv Survev CohaborationI 
120051) , the Panoramic Survey Te lescope and Rapi d 
Response System^^ (PanSTARRS; IKaiser et al.l[200l . 
the P hysics of the Accelerating Universe survey^* 
(PAU; iBenftez et all 120091), the Lar ge Synoptic Sur- 
vey Telescope^^ (LS ST; iTvsonI l200l . and EUCLID^" 
(|Beaulieu et al.ll20Tot) . 

A number of previous works have analyzed and re- 
ported the cosmological constraints from the galaxy clus- 
tering of the imag in g surveys (e . g.. iTb g piark ct al. 20021 
Blake etal.1 I2007t iRoss et all 120081; ISawangwit et alj 



20091; iThomas et al.l I2010L 120111 : ICrocce et al.l l20TlF 
Ross et al.l 120111 )^ but there have been only a 



few published works on the BAO measu rement 
(IPadmanabhan et all 120071: ICarnero et all 1201 ID . Our 
goal in this paper is to design a robust method for mea- 
suring the location of BAO in the angular power spec- 
trum of imaging surveys and apply it to the final imag- 
ing data set of t he Sloan Digital Sky Survey HI (SDSS; 
I York et al.ll2000l) . 

We use the DR8 imaging catalog of SDSS-III that in- 
cludes photometric redshifts of luminous galaxies (here- 
after, 'LGs') betw een 0.45 <z< 0.65 over - 10, OOOdeg^ 
()Ross et al.l 1201 ID ; the spectroscopy from the SDSS- 
III Baryon Oscillatio n Spectroscopic Survey (BOSS; 
lEisenstein et al.ll2011|) is used to create a training sample 
and therefore to estimate the true redshift distribution 
of photometric galaxies. The angular power spectra are 
generated from the data using an optimal quadratic esti- 
mato r, as presented in detail in one of our companion pa- 
pers, [Hol^i^l] ()201lD . Utilizing the estimated true red- 
shift distribution, we construct a theoretical BAO tem- 
plate for the angular power spectrum, fit to the loca- 
tion of the observed BAO feature in the angular power 
spectrum, and derive the angular diameter distance to 
z = 0.54 in a manner independent of underlying dark 
energy models. 

The angular clustering of galaxies contains more cos- 
mological information than the scale of BAO. Redshift 
distortions on very large scales and the overall shape of 
the power spectrum (e.g., the matter-radiation equality 
feature) can provide additional information. However, in 
this paper, we take a very conservative approach and use 
only the most robust probe, the location of B AO, while 
excluding most of the non-BAO information. IHo et all 
(j2011l Paper I) presents a more extensive study; it in- 
cludes information of the full power spectrum and derive 
cosmo logical constraints. In p arallel, another companion 
paper, Ide Putter erall (|201ll) . measures the mass bound 
on the sum of neutrino masses using the same power 
spectra. 

This paper is organized as follows. In § [21 wc briefly 
summarize the imaging data. In § [31 we summarize the 
method used in Paper I to generate angular power spec- 

26 www.darkcncrgysurvey.com 

27 www.pan-starrs.ifa.hawaii.edu 
2* www.icc.cat/pau 

23 www.lsst.org 

30 www.sci.esa.int/cuclid 
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Fig. 1. — The true rcdshift distribution estimated for our pho- 
tometric redshift galaxies for the four redshift bins: CMASSl for 
0.45 < Zph < 0.5, CMASS2 for 0.5 < z^^ < 0.55, CMASS3 for 
0.55 < 2ph < 0.6, and CMASS4 for 0.6 < Zph < 0.65. The median 
and the mean of the combined galaxy distribution is 0.541 and 
0.544, respectively. 



tra. In §111 we describe details of the BAO fitting method 
used in this paper. In § [51 we test our method and as- 
sumptions using mock data. In § [6l we present an analy- 
sis of the DR8 imaging data and report the best fit angu- 
lar diameter distance to sound horizon ratio at z = 0.54. 
We show the robustness of our result and show the effect 
of correcting for the observational systematics. In § |7l 
we discuss constraints on various cosmological parame- 
ters when our BAO measurement is combined with the 
WMAP7 data and other BAO measurements. Finally, in 
§ [51 we summarize the results in this paper. 

2. DATA 

We use th e imaging data of th e eighth and fi nal data re- 
lease (DR8; lAihara et al.llMll) of SDSS-III (lYork et al.l 
|2000[) that is obtained by wide-fiel d CCD phot ometry in 
five passbands (u g, r i, z) (see lFukugita~et al. 199^; 
IGunn et al.|[l998l [200l [Pier et al.ll2003l for more tech- 
nical and data realease details). We use the p hotometric 
redshi ft catalog constructed as described in iRoss et all 
()201H) '^^ This catalog is selected fro m DR8 using the 
same criteria as the SDSS-III BOSS (jEisenstein et al.l 
120 111) targets selected to have approxima tely constant 
stellar mass (CMASS: iWhite et al.ll2011| ). Photomet- 
ric redshifts and probabilities that an object is a galaxy 
were obtained using a training sample of 112,778 BOSS 
CMASS spectra [to be released with Data Release 9 in 
July 2012] 32. The final catalog covers 9,913 deg^ of sky 
and consists of 872, 921 galaxies between 0.45 < z < 0.65, 
which is an improvement in the s urvey area compared 
to the MegaZ-LRG DR7 catalog (|Thomas et al.l (20111 
723,556 objects over 7, 746 deg^ for 0.45 < z < 0.65). 

The estimated photometric error, CTz^j^, increases from 
0.04 to 0.06 ove r the redshift range (See Figure 10 of 
iRoss et all 1201 1[ ). We define four photometric redshift 
bins, with widths similar to the photometric error, re- 



TABLE 1 

The four photometric redshift bins. 





bins 


Zph range 






-^median 


^mcan 




CMASSl 


0.45-0.50 


•2u,m 


0.043 


0.474 


0.475 




CMASS2 


0.50-0.55 


258,736 


0.044 


0.523 


0.526 




CMASS3 


0.55-0.60 


248,895 


0.052 


0.568 


0.572 




CMASS4 


0.60-0.65 


150,319 


0.063 


0.617 


0.621 




Total 


0.45-0.65 


872,921 




0.541 


0.544 



Note. — A'gai is the effective number of galaxies 
after weighting each object by Psg, which is the proba- 
bility that an object is a galaxy. The value of cr^pj^ is 
the dispersion in redshift for each photo-z bin. 



ferred to as CMASSl, CMASS2, CMASS3, CMASS4. 
Table [H and Figure [1] show the distribution of the effec- 
tive^'^ number of galaxies in each redshift bins. Due to 
the extensive training sample, our determination of the 
redshift distribution is expected to be quite accurate; for 
example, based on the Jack-knife resampling of the train- 
ing sample, we estimate the error on the mean/median 
of the distribution of each redshift bin to be less than 
0.5%. The median and mean of the combined galaxy 
distribution arc 0.541 and 0.544, respectively. 

3. OPTIMAL QUADRATIC ESTIMATOR OF 
ANGULAR POWER SPECTRA 

The auto and cross angular power spectra of the 
four redshift bins were generated using the optimal 
quadratic estimator in Paper I, to which we refer the 
readers for more details of the optimal quadratic csti- 
mator (OQE, also se e Scliak 1998; Tegmark et al. 199J; 
iPadmanabhan et al.l [20031 12007). To summarize, we 
parametrize the power spectrum with 35 step-function 
band powers and write the data covariance matrix as 



(/3) 



(1) 



where Si is the galaxy overdensity at the i pixel at 
a given redshift bin, is the band power for a wave 
number bin /3, C^^^ is the derivative of C with respect to 
the band power pp, and N is the shot noise contribution 
to the covariance matrix, while Sj is the Kronecker delta 
function. Assuming Si is Gaussian-distributed, requiring 
the estimator to be unbiased and to have a minimum 
variance, we derive a band power estimator 



^[3-1 



StQ-iCh)C-iS -b. 



where b-y is the 
itrp-iC^-^'C-iiV]), i.e. 
the Fisher information matrix 
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contribution from 
the shot noise, 



N 
and 



F, 



1 



Pi 



The variance of the band power is derived by 

Cov[pi3,p^] ^ F, 



Pi 



(2) 

(i.e., 
F is 



(3) 



(4) 



Available at http : //portal . nersc . gov/proj ect /boss/ galaxy /photozp -y^e weight each object by the probability that an object is a 
We use the redshifts available through MJD 55510 galaxy. 
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Fig. 2. — The measured angular power spectra for the four red- 
shift bins. The solid lines show the best fits derived in § [6] The 
dashed line is the template for CMASS2 which is rescaled for clar- 
ity, with redshift distortions assuming a galaxy bias of 2. The 
dotted line is the template without redshift distortions. 



The expected value of band power, {pjs), is related to 
the power spectrum p{£') a t an integer £' by the band 
window function Wpe' (e.g.. IKno3ll999f ): 



dp{e) 



p{£') = W0tp{£'), 



and 



(5) 



(6) 



where C*-^ ^ is the derivative of C with respect to the 
power at an integer wavenumber £' . In this paper, we 
use the band power pp and the covariance matrix of the 
band power derived in Paper I. The auto power 

spectra are shown in Figure [2j The cross-power esti- 
mates between different redshift bins were generated in 
a similar manner and used only for determining covari- 
ance between different redshift bins. 

4. METHODS 

4.1. Outlines of the fitting method 

Our goal is to robustly measure the location of the 
BAO scale, and therefore the distance scale Da(z), while 
minimizing possible effects from the assumptions made 
on the nonlinear or/and observational effects as well as 
cosmology during the fitting procedure. For example, 
it is non-trivial to properly model the evolution of the 
broadband power on the nonlinear scales, whether it is 
due to structure growth, redshift distortions, or galaxy 
bias. Therefore, we aim to exclude as much non-BAO 
signal as possible by marginalizing over the effect of the 
smooth broad-band power. We measure the location of 
the BAO feature by fitting the observed auto (band) 
power spectrum Cobs,zi {(-Y^ with the following fitting for- 
mula using a template power spectrum Cm.7.-X^/oi)'- 

Cobs,z. W = B..(£)a„.z. Wa) + (7) 

Note that we switched the notation from to C(t). 



where a, Bzi{£), and Az-{i) are fitting parameters. The 
functional form for . and A^. is discussed in § 14.21 The 
parameter a measures the angular location of the BAO 
relative to that of the fiducial cosmology. That is. 



Lhs/hd = [DA{z)/rs]ohs/[DA{z)/r 



j,(r(fc,z))-/3j, (r(fc,z)) p] 



(8) 

where [DA{z)/rs]fid is the fiducial angular location of the 
BAO in the template and [DA{z)/rs]ohs is the measured 
angular location of the BAO. A value of a > 1 sug- 
gests that the observed angular location of the BAO is 
smaller than that of the fiducial cosmology. For each red- 
shift bin, Zi, free parameters B^^t} and Az-{i) account 
for the smooth modification of the power spectrum from 
the template due to nonlinear structure growth and any 
scale-dependent bias. Finally, we use power spectra for 
CMASSl, CMASS2, CMASS3, and CMASS4 simuhane- 
ously and fit for a universal a while marginalizing over 
and A^- independently at each redshift bin. 
The template Cm,zi(^/a) is constructed from the 
2-dimensiona.l projection of 3D power spectrum 
(jFisher et al.l 119941: (Padmanabhan et al.ll2007D . Includ- 
ing linear redshift distortions, 

C,„,z.W dkeP^{k,z,) 

dN , D{z) 

dz biz) — 

dz ^ 'd{z = 0) 

where r(fc, z) = k{l + z) Da, fid{z), dNzJdz is the normal- 
ized, true redshift distribution of galaxies for the corre- 
sponding i*'* photometric redshift bin (Figure [J), is 

the spherical Bessel function, is the second derivative 
of the spherical Bcssel function with respect to r(fc, z), 
6(z) is a fiducial galaxy bias, /3 is the fiducial redshift dis- 
tortion parameter'^^, and D{z) is the linear growth rate. 
Due to the projection, redshift distortions significantly 
affect the broad-band shape of the power spectrum only 
for I < 30; s ee the dashed and dotted lines in Figured] 
(also see e.g.. lNock et~aIll2010D . 

Knowing dN^-fdz precisely is critical for constructing 
the correct BAO location in the template. Note that, 
thanks to the extensive spectroscopic training set from 
the BOSS CMASS galaxies (112,778 objects), we have 
an excellent determination of the redshift probability dis- 
tribution for our photometric luminous galaxy samples. 
Once we calculate the template Cm,zi(0 for the fiducial 
cosmology, we rescale the wavenumber with a to gener- 
ate Cm.zi(^/a) that fits the observed power spectrum in 
Equation [7] while marginalizing over the free parameters 
Bz,ii) and Az^{e). 

The term P,„ is generated by degrading the 
BAO portion of the fiducial, linear power spec- 
trum with a nonlinear damping parameter — 
7.527[D{z)/D{0)]h-^ Mpc to mimic the nonlinear evo- 
lution of the jBAO due to the structure growth 
(jEisenstein et al.ll2007| ): 

P,n{k, Z,) = [Piinik) - P„w(fc)] CXp [-k''^,n{zf /2] 

+ Pnw(fc), (10) 

where Pun is the linear power sp ectrum at z = derived 
from CAMB (jLewis et al.ll2000[) and Pnw is the nowiggle 
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form, i.e., th e power without B AO, calculated using the 
equations in lEisenstein fc Hul ()1998f) . The smoothing of 
the BAO is dominated by the width of the underlying 
redshift distribution, and the exact choice of Em(z) does 
not have a significant impact. For clarity we refer to 
Cm,zi(^), instead of Pm, as a 'template' (angular) power 
spectrum and Pm as a 'base' power spectrum in this pa- 
per. We use a fiducial cosmolo gy similar to the WMAP7 
results (jKomatsu et al.l [2011[ ) to define Pm, -DAfid(-z), 
and etc: = 0.274, = 0.726, h = 0.7, = 0.046, 
Us = 0.95, and erg = 0.8. These values produce the 
fiducial sound horizon scale, = 153.14 Mpc, based on 
lEisenstein fc Hul ()1998D , and the fiducial angular location 
of the BAO, [DAiz)/rs]fid = 8.585 at z = 0.54.3^ 

Note that, in equation |9l we are assuming that rescal- 
ing the sound horizon and angular diameter distance is 
equivalent to rescaling i (i.e., the second equality of Eq. 
[8]), which we call 'a model'. This will be a reasonable 
approximation when the thickness of the redshift distri- 
bution that is projected on the 2D celestial surface is 
much larger than the scale of the clustering, i.e. in the 
limit of the Limber approximation. Based on the mock 
test in §[5] and the robustness test in §|6l the a model ap- 
pears to be a good enough approximation for our choice 
of fitting range and parameterization. 

The measured band power Cobs,zi (^) has a contribution 
from a range of wave numbers, which is described by a 
window function (Eq. [B]). Due to the large and contigu- 
ous survey area, the window function is sharply peaked 
with little correlation with neighboring bands (Figure [3]) • 
We account for this window function effect in the fitting. 
That is, for each wavenumber band Ia for the redshift 
bin Zi, 

C'obs,Zi(^A) = 
(.600 

X [B,^{£')C^^,^{e' /a) + A,^ie')]w,i£'\£A)de\n) 



4.2. A choice of B{1) and A{i) 



0.3 



0.2 



0.1 




Fig. 3. — Window function of CMASS2, as an example. 

Using the exact calculation of dr ag epoch rather than the 
fitting formula in lEisenstein fc Hul {1991), we find rs = 149.18 Mpc 
for our fiducial cosmology, therefore lDj^{z)/rs]fid = 8.813 at z = 
0.54. 



There are a few considerations to make when choosing 
the optimal parametrization for B{£) and A{£). We want 
B{£) and A{£) to be flexible enough to model and remove 
the broad-band shape of the power spectrum. With the 
flexibility in B{£) and A{£), we also gain some toler- 
ance on a possible, small difference between the BAO 
feature in the fiducial template and the observed fea- 
ture by trading power between i?(^)Cm,zi (^) s-i^d ^(^)- 
On the other hand, an arbitrarily flexible B{£) and A{£) 
will undesirably mimic BAO even with the no-BAO tem- 
plate. An extensive test of the parametrization for B{£) 
a nd A(£) fo r the template fltting method is discussed 
in lSeo et al.l (I2008D . where they allow a large number of 
free parameters for B and A based on the spherically- 
averaged power spectra from A'^body realizations which 
correspond to a spectroscopic survey. With our photo- 
metric redshift uncertainty, a projection of the BAO from 
different distances introduces an additional damping in 
the feature, leaving higher harmonics other than the flrst 
much less distinct relative to the noise level. Therefore 
we are forced to limit the flexibility i n B(£) an d A(£) 
more strictly than the case assumed in lSeo et al.l ()2008l ) 
while still making sure that the BAO scale is correctly 
recovered. We choose a revised fltting range so that the 
broadband is well modeled despite the smaller number of 
B{£) and A(£). 

Based on tests with mock catalogs (§[5]), we use a fltting 
range 30 < / < 300 and a linear function in £ for B^. and 
a constant A^. (abbreviated with 'AOBl' hereafter): 



A, 



B, 



(12) 
(13) 



Therefore, for four redshift bins, we flt for a total of 13 
parameters including a. 

4.3. Parametrizing Da{z) 

The observed location of the BAO peaks in a power 
spectrum is determined by the angular diameter distance 
at each redshift. In a sample selected using photomet- 
ric redshift, the location of the BAO feature for a given 
redshift slice depends on an integration over the broad 
range of the true redshift distribution, rather than a dis- 
tance at a single redshift, as evident in equation [HI In an 
ideal case, we may attempt to constrain the redshift de- 
pendence of Da{z) over the entire s pectroscopi c redshift 
range in a non-parametric way (e.g.. lPercival et al., 2010[) 
in the fltting process. In reality, the Fisher matrix anal- 
ysis predicts an error of ~ 4% on a single measurement 
of Da{z) for all our four redshift bins combined. Also 
while the expected true redshift distribution has long 
tails (0 < < 1), the distribution tends to peak sharply 
for each photometric bin so that most of the informa- 
tion for a given bin is well concentrated within icr^pj^ . 
Third, for a reasonable range of cosmology, the shape of 
Da{z) does not evolve significantly between z = 0.45 and 
z = 0.65. We therefore expect little gain for designing 
our analysis to measure multiple Da{z)s, i.e., the evolu- 
tion of Da{z) as a function of redshift. We instead design 
our fitting method to measure a single, more precise dis- 
tance measurement at the redshift that contributes the 
most information. Based on the median and mean red- 
shift of the photometric sample, we assign our measure- 
ment of Da to z = 0.54. 
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In detail, wc assume DA.ndiz), given by the fiducial 
eosmology, and set 



Da{z) = aDAfid{z)- 



(14) 



That is, we fix the shape of the Da{z) to be the same as 
DA.fid{z) and measure the amplitude of Da{z). 

4.4. Clustering evolution of lumious galaxies 

In generating the template Cin,zi(^), we need to make 
a prior assumption on the evolution of the galaxy bias of 
LGs and the linear growth rate (Eq. [9]). We consider two 
extreme cases of the galaxy clustering evolution: first, we 
assume that the overall clustering, b^D^, does not change 
with redshift, which we call as ^ con-cluster\ Second, we 
assume that the bias does not change with redshift, which 
we call as ^con-bias' . The two cases make little difference 
in the final best fit of a, mainly because the expected true 
redshift distribution sharply peaks within ztcr^pj^, com- 
pared to the galaxy clustering evolution. Note that, by 
marginalizing over B^. at each photometric redshift bin 
Zi, we take into account the evolution of galaxy clustering 
across different redshift bins whether we use 'con-cluster' 
and 'con-bias'. As a default, we fix & = 2 inside Cm,zi(^) 
(i.e., ^ con-bias\ and therefore the best fit B^^ can be ap- 
proximately interpreted as b'^[zi). 

5. TESTING THE METHOD 

Before applying our fitting method to the real data, we 
want to validate, using mock catalogs, that our fitting 
method returns an accurate estimate of the BAO scale. 
In other words, we want to check that neither our process 
of deriving optimal quadratic estimators of band powers 
nor our fitting method biases the measured BAO scale. 

5.1. N-body mocks 

As explained in Paper I in detail, we generate mock 
catalogs of our imaging data making use of th e 20 
CMASS mocks constructed bv IWhite et~all ()201lD . We 
call these 'A^-body photoz-mocks'. The cosmology used 
for generating these mocks is the same as our fiducial cos- 
mology. The comoving volume of the original CMASS 
mock is [l.bh^^ Gpc]'^ and, to build A'-body photoz- 
mocks, we extract an octant of a spherical shell between 
r = 1.33/1-1 Gpc(z = 0.5) and lA5h~^ Gpc(z = 0.55) 
from the origin (one corner of a simulation box) and 
project galaxies along the radial direction without in- 
troducing photometric redshift errors. For simplicity, we 
do not include redshift distortions, which will be visi- 
ble only on very large scales (Figured]), nor the effect 
of the mask in generating these mocks. The resulting 
power spectrum of the projected field has a BAO feature 
that is quite similar to that expected for CMASS2(i.e., 
0.5 < z < 0.55). Each of the resulting A-body photoz- 
mocks spans a 7r/2 rad^ (= 5157 deg^) and contains 
~ 125,000 galaxies. While the number of galaxies is 
smaller than CMASS2, the amplitude of power spectrum 
is boosted overall by almost the same factor due to the 
thinner redshift slice: the signal-to-noise ratio per mode 
of each mock is therefore similar to CMASS2. (see Fig- 
ure H]) . Therefore the photoz-mocks serve as reasonable 
mocks for the observed power spectrum. 

We repeat the procedure by placing an origin at eight 
different corners of each simulation box and generate 
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Fig. 4. — The red circles with error bars show a power spectrum 
averaged over 20 Af-body mocks for the same line of sight. The 
black points show the power spectrum of CMASS2. The dotted 
lines in the top panel show the shot noise contribution in both 
cases. The solid line in the bottom panel show the expected BAO 
feature. 
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photoz-mocks for different lines of sight. For each line of sight, the 
power spectrum is averaged over 20 A^-body photoz-mocks. Note 
that the different lines of sight are correlated. 



eight sets (i.e., eight lines of sight) of 20 imaging mocks, 
i.e., a total of 160 mocks. Note that the eight lines of 
sight from each simulation box share a portion of volume 
and therefore are not independent of each other. We gen- 
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Fig. 6. — The top panel shows one of the Gaussian CMASS mocks for CMASS2 (magenta circles) in comparison to the real data (black 
squares). The solid lines show the best fit model we used to generate the mocks. The other three panels show the distribution of the 500 
best fit a values of the Gaussian CMASS mocks. Top-right panel: the template is generated using the same cosmology as one used for 
generating the mocks (i.e., fim/i^ = 0.134). The histogram shows the result using AOBl over 30 < I < 300. The black point with error bar 
shows the mean a — 1 and the error associated with the mean. The black dot-dashed line near the center shows the expected best fit a — 1 
for each cosmology. The red and magenta (only in the top second panel) dotted lines show the average error derived from the 500 mock 
values for the range for Ax^ = ±1 and the range for the 68.3% of the likelihood, respectively. The blue dotted line shows the range that 
contains 68.3% of the distribution of the 500 best fit as. The three values arc quite similar. Bottom-left: the distribution of the best fit a 
of CMASS mocks when the template is built using Qmh^ = 0.148. Bottom-right: using the template built assuming Q,nh^ = 0.127. 



eratc the template power spectrum based on the galaxy 
distribution of the mock catalogs using Equation [9l 

To better detect a possible bias on a when using our 
fitting method, we increase the signal-to-noise ratio by 
averaging many power spectra. We average power spec- 
tra of the 20 mocks for each configuration (i.e., each line 
of sight) and fit for a. For the eight configurations, us- 
ing 'AOBl' (i.e., with Bq + Bi£ and Aq), we show the 
derived best fits and the 68.3% range of the likelihood 
in Figure [S] The distribution of eight a values appears 
slightly wider than what one would expect from a Gaus- 
sian distribution. However, it is not appropriate to com- 
pare this result with a Gaussian case, as the number of 
sets (i.e., eight) is too small to estimate the underlying 
distribution and the eight configurations are not inde- 
pendent from each other but share some portion of their 
volumes. Overall, we do not see an obvious indication ei- 
ther that the best fit a from our OQE estimator is biased 
or that our fiducial fitting method introduces a bias. 



5.2. Gaussian CMASS mocks 

We next test the accuracy of our fitting method using 
sets of mocks that have exactly the same noise prop- 
erties as the real data. We generate many Gaussian 
mock power spectra using the covariance matrix for the 
real data described in § 16.11 In detail, we find the best 
fit band powers for the power spectra measurements of 
CMASSl, CMASS2, CMASS3, and CMASS4 over the 
entire wavenumber range (i.e., 1 < / < 600), while delib- 
erately using a slightly different choice of B{1) and A{1) 
(A2B0) than our fiducial method; a is fixed to be unity 
during this process. This approach allows us to con- 
struct theoretical power spectra for the fiducial cosmol- 
ogy while taking into account the realistic broad-band 
shape of LGs. Using the four best fit band power spec- 
tra and using the covariance matrix for the real data, 
we generate 500 sets of Gaussian CMASS mocks. For 
each set, the four CMASS mock power spectra therefore 
mimic CMASSl, CMASS2, CMASS3, and CMASS4 not 
only in terms of clustering but also in terms of the covari- 
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ance among them. One of the mock CMASS2 is shown 
in the top left of Figure |6] in comparison to the real data. 

We then apply our fiducial fitting method to the 500 
CMASS mocks. Since the mock data and the template 
use the same fiducial cosmology, we expect the average 
value of a, if unbiased, to be unity. The top right panel 
of Figure [5] shows the pdf distribution of the 500 best fit 
a values when we fit to the mock data over 30 < ^ < 300 
using AOBl. The mean and the standard deviation of 
a — 1 is (0.10 ± 6.60)%. After rescaling the standard 
deviation by the square root of the number of samples, 
the mean and the error associated with the mean value 
of a is (0.10 ± 0.30)%, i.e., unbiased within 1 — a. De- 
creasing the number of free parameters to AOBO causes 
a biased estimation of a as large as 2% in a. We find 
AlBO also gives an unbiased result, while we find that 
paramctrizations with more free parameters result in a 
non-negligible degree of bias on the BAO location. 

The distribution has significant 'tails' as shown in Fig- 
ure [6] As a result, the standard derivation (6.60%) is 
much larger than the range that contains 68.3% of the 
distribution: a - 1(%) = O.lttl- That is, the distri- 
bution of a appears wider than the Gaussian one. We 
have excluded a few catastrophic outliers (~ 0.8% of the 
samples for the fiducial case) that show a — 1(%) > 40 in 
deriving these statistics. Obviously, these cases are not 
fitting to the BAO feature. Some of these catastrophic 
outliers show a factor of a few larger errors on a than 
the rest of the samples, while others show a reasonable 
error associated with the best fit. We find that the re- 
duced of these catastrophic outliers are not necessarily 
large. While we exclude these extreme outliers so that 
the statistics of the distribution are not dominated by 
these occasions, we note that, in analyzing a real data, 
we expect to derive a very wrong value of BAO scale 
more likely than would be expected were this distribu- 
tion perfectly Gaussian. 

For the real data, we have only one realization and 
therefore need to check that the error we quote for the 
data will be closely approximating the 68.3% range of the 
sample distribution if we had more than one sample. For 
each CMASS mock set, we derive errors associated with 
the best fit a by applying Ax^ = ±1 or by deriving the 
width that contains 68.3% of the likelihood. On average, 
the resulting error on a for an individual CMASS mock is 
a-l(%) = O.lttI and 0.lt^:7 for Ax^ of±l (dotted red 
lines in Figure |6]) and the 68.3% width of the likelihood 
(dotted magenta lines) , respectively; the latter is slightly 
larger than the former. These are reasonably similar to 
the 68.3% range of the best fit distribution of the mocks, 
which was a - 1(%) = O.lttl (dotted blue lines). We 
therefore will quote the 68.3% range of the likelihood 
surface as our formal error for the real data. 

5.3. Variations in the template 

In Section 15.21 we assumed that the fiducial cosmol- 
ogy used for constructing the template matched the true 
cosmology, and we could recover an unbiased result in 
this situation. In fact we also need to confirm that our 
method is unbiased when the cosmology used in the anal- 
ysis does not match the true cosmology. The shape of 
the BAO feature is determined by the matter density and 
the baryon density, which is well measured by the cur- 



rent CMB observations. iSeo et ahl ()2008D investigated 
how the results of the template fitting depends on small 
deviations in Pm and found that the effect is negligible 
(less than 0.02% bias on a). As explained in § 14.21 our 
method in this paper uses a smaller set of free parame- 
ters due to the lower signal-to-noise level of the real data 
and, as a result, it is possible that the fitting result is 
more sensitive to the deviation in P^- We therefore re- 
visit this issue. Since n^/i^ is better measured than fim^^ 
by the CMB data, we only vary flmh'^ from the fidu- 
cial value, while th e current CMB const raint on flmh^ 
is 0.1326 ± 0.0063 (jKomatsu et al.llMl . In order to 
leave the shape of the angular diameter distance to red- 
shift relation unchanged, we hold 0,,^ fixed and vary h 
accordingly. 

The bottom left panel of Figure IH] shows the distribu- 
tion of the best fit a of 500 Gaussian CMASS mocks 
when the template is built using ftmh^ = 0.148 (i.e., 
10% away from ilmh? = 0.134 used for the mocks). We 
again use the fiducial, AOBl parameter set and a range 
of 30 < Z < 300. We find that using a smaller number 
of parameters, i.e., AOBO, makes the result more vul- 
nerable to the variations in the template. Based on the 
sound horizon scale and h in this cosmology, we expect 
a — 1(%) = 1.6 (the dashed vertical line). The average of 
the best fit is a — 1(%) = 1.96±0.31, and therefore we re- 
cover the correct BAO scale. The bottom right panel of 
Figure ini shows the distribution of a using ilmh? = 0.127, 
i.e., —5% away from the true value used for generat- 
ing the mocks. We expect a — 1(%) = —2.10 while we 
measure —1.84 ± 0.27; we recover the expected value. 
We also have tested AlBO: we find a moderate bias for 
Vlmh? = 0.148 due to asymmetric tails, while the bias 
is overall less than 0.6%. This parametrization therefore 
would be a reasonable choice as well given the level of 
signal to noise of our data. 

In the next section, we will apply the same test to 
the real data and show that the measured BAO scale 
does not change as a function of the fiducial cosmology 
assumed for the template. In addition to the assumption 
wc made for P,„ , we have also assumed a fiducial relation 
of angular diameter distance to redshift. We will show 
that we recover the same result over a range of angular 
diameter distance to redshift relationship. 

Summarizing our mock tests, we find that the angu- 
lar p ower spectra produced by our OQE code (|Ho et al.l 
1201 It ) show no strong sign of a bias in its BAO feature 
and that 'AOBl' is a good choice of parametrization over 
30 < / < 300 for our data quality and therefore returns 
the correct BAO scale within 0.3% over a reasonable 
range of variations in our assumption. 

6. RESULTS: DR8 IMAGING DATA 

6.1. Building the covariance matrix 

We use Gaussian covariance matrices calculated for the 
auto power spectra of the four DR8 redshift bins using 
equation m As shown in Paper I, we do not find an ob- 
vious indication that the Gaussian covariance matrix for 
the OQE estimator underestimates the true error of the 
2-dimensional projection of the nonlinear galaxy field. 

ISeo et all 1200^) uses an Af-body volume of 320h-'' Gpc^, as 
a comparison. 
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Also lTakahashi et alj (|2011[ ) and lNgan et"all (1201 1[ ) have 
shown a neghgible effect of non- Gaussian errors on the 
BAO measurement in muUi-parameter fitting. We there- 
fore retain the Gaussian field assumption when deriving 
the covariance between different redshift bins as well: 
we use the cross-power spectra between different redshift 
bins, while taking the window function due to the survey 
mask into account (see Paper I for more details). 

6.2. Best fit angular location of BAO 

We apply our fitting method to the DR8 imaging data 
and constrain the angular location of the BAO. Figure 
[7] shows our best fit result using the combinations of 
CMASSl, CMASS2, CMASS3, and CMASS4. The red 
lines/points show the best fit C; with the BAO template 
in comparison to the measured data (black squares with 
error bars). The range of red points show the range of 
the fitting, i.e., 30 < Z < 300. In the figure, we denote the 
best universal fit a with the associated errors that corre- 
spond to the 68.3% range of the likelihood distribution: 
we derive a - 1(%) = 6.61^4:^^. The reduced at the 
best fit is 1.20 for 87 degrees of freedom, and the proba- 
bility of having a reduced x value that exceeds this value 
is 10%. The left panel of Figure [8] shows the resulting 
surface along a when marginalized over other parameters 
(red line). Note that, due to the oscillatory feature of the 
BAO both in the data and the template, there are local 
minima around the global minimum of x^- As implied 
in the figure by the extent of the red line, when we de- 
rive the 68.3% range of the likelihood, wc only include x^ 
over 0.08 < a < 1.4, avoiding the local minima beyond 
this range. The right panel of Figure [5] shows a stacked 
Ci/Ci,sm of the four panels of Figure [71 We perform this 
stacking procedure as follows. To better visualize the 
BAO feature we measured, we shift the wavenumbers of 
the four power spectra by DAizmodmn) / Da{z = 0.54), 
where DA(-Zmedian) is the median redshift for each red- 
shift bin. We combine the four band powers, re-bin 
the combined data while inversely weighting each band 
power by its error. The solid red line is the best fit for 
CMASS2 after its wavenumber is rescaled to mimic a 
result at z = 0.54. 

Interpreting a — 1(%) = 6.61 requires our determi- 
nation of the redshift to which this measurement corre- 
sponds. Strictly speaking, the best fit value of a rep- 
resents a constant ratio of the observed DA{z)/rs to 
the fiducial [D a{z) / r s]fLd assumed in the template. The 
black solid and dashed lines (with a shade) in Figure [HI 
show what the best fit and the 1 — a error on a imply in 
this strict interpretation. However, although the redshift 
dependence of Da{z) we assume spans z — 1, most 
of the galaxies are within z = 0.45 and 0.65 with a peak 
of the distribution near 0.5 < z < 0.55. Therefore, it 
is reasonable to consider that the best fit a[D a{z) / r s]fid 
represents DA{z)/rs near z = 0.5—0.55. The median and 
the mean of the weighted galaxy distribution are 0.541 
and 0.544, respectively. We therefore adopt z = 0.54 as 
the characteristic redshift that our BAO measured scale 
represents. 

To show that the best fit BAO scale indeed does not 
depend on the cosmology we assume for [DA{z)/rs]fid, 
we repeat our fitting with templates constructed using 
different cosmologies. In detail, we vary the equation of 
state of dark energy, w, by ±0.2, such that [DA{z)/rs]tid 
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Fig. 7. — The best fit result using tiie combinations of CMASSl, 
CMASS2, CMASS3, and CMASS4 and using AOBl. We derive 
a — 1 = 6.609 j|4'gg%: i.e., 4% deviation from the fiducial value 
based on WMAP7. The black data points in the four panels 
with error bars show the measured C'l divided by a smooth fit at 
CMASSl, CMASS2, CMASS3, and CMASS4. The red Hues show 
the resulting best fit C; and the red circles show the best fit band 
power Ci after the window function effect is considered. 
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Fig. 8. — Left: the surfaces along a for Figure [7] when marginalized over other parameters (red line). The vertical dotted lines show 
the best fit a and the 1 — a range. The blue line shows the surface when the BAO is removed from the template. Right: a stacked 
Ci/Ci of the four panels of Figure [7] To bettor visualize the BAO feature we measured, we shift the wavenumbers of the four power 
spectra by Da (^median ) A = 0.54), re-bin the combined band powers while inversely weighting by errors. The solid red line is the best 
fit for CMASS2, after its wavenumber is rescaled to ^ = 0.54. 



at z ~ 0.54 varies by ^ 3.4 — 3.7%, i.e., slightly less than 
the 1 — cr range associated with the best fit a. In other 
words, we are testing the consistency of our answer by 
varying the template by 1 — cr from the fiducial case 
of a. Using a template with w = —1.2 , the best fit 
gives a - 1 = 2.98ttil% and, with w = -0.8, a - 1 = 
10.2115-02%, 

The top panel of Figure [H] shows the best fit DA{z)/rg 
(= aDA{z) /rg^fid) using the three different templates. 
The solid lines show the best fit and the dotted lines 
show a 1(7 range of DA{z)/rs- The bottom panel displays 
ratios of the best fit DA{z)/rs using no-ACDM templates 
with respect to the best fit DA{z)/rs using our fiducial 
ACDM template. From the top and bottom panel, one 
sees that the three different templates have a very similar 
shape in Da{z) over z = 0.45 — 0.65, once the absolute 
difference is absorbed into a. The three templates re- 
turn virtually the same DA{z)/rs at z = 0.54: they are 
consistent within 0.3%. We therefore quote the best fit 
using our fiducial ACDM template as our official mea- 
surement: DA{z)/rs = 9.212toill at z = 0.54. Using 
the current WMAP7 constraint on the sound horizon at 
drag epoch, 153.2 ± 1.7 Mpc, we derive angular diameter 
distance Z?4(z) = 1411 ± 65 Mpc at z = 0.54 ^s. Table 
[2] summarizes our best fit BAO location and the derived 
distance scale. 

We further test the robustness of our result by con- 
structing templates using various cosmologies. Figure 
[TOl shows the best fit DA{z)/rs at z = 0.54 assuming 
wCDM, oCDM, ACDM, and assuming different values 
of ^mh'^- For the range of cosmologies we have investi- 
gated in this paper, the best fit varies less than 1% in the 
acoustic scale while the la error is ^ 4.7%. The errors 
vary slightly more than the variations in the best fit, es- 
pecially when the template cosmology deviates substan- 
tially from the concordance cosmology. We also test a 

^® U sing the exact integrat ion rather than the fitting for- 
mula in lEisenstein fc Hul llT998tl . = 149.18 Mpc and Da{z) = 
9.456lo 4i5 at z = 0.54. 



different parametrization than the fiducial choice, AlBO, 
which has marginally passed the mock test (i.e., a likely 
bias of ~ +0.6% on a based on the result in i^ l5.3p . AlBO 
gives a — 1 = 7.5(%), which is consistent with the fiducial 
result within 1%. 

If we remove the BAO in the template (the blue line in 
Figure E]), we essentially fail to constrain a. This means 
that the flexibility in our fitting is sufficient that the 
broadband shape information cannot constrain a. There- 
fore, we conclude that our measurement of DA{z)/rs is 
mainly from the BAO information. 

The p recision of our measu rement is much better than 
that of iCarnero et al.l (j2011| ) , where they measure the 
BAO location within 9.7%. The discrepancy arises partly 
because of the larger survey-area coverage in this work, 
but the m ain difference i s due to the difference in the 
modeling. ICarnero et al.l ()2011[ ) assume a quite general 
model without utilizing the true redshift distribution of 
their photometric galaxies, and account for the various 
scenarios of the deviation between their model and the 
real data in their error budget. Our method is, on the 
other hand, a template-based approach utilizing the true 
redshift distribution of our photometric galaxies that is 
quite well determined by the extensive training set and 
generate a template with a precise BAO location given 
the redshift distribution and cosmology. Our method 
therefore is quite immune to the most of the systematic 
errors they account for. We observe approximately 1% of 
variations depending on the choice of fiducial cosmology 
and parametrization, which has very small effect when 
added quadratically to our current error of 4.7%. 

6.2.1. Significance of detection 

A reasonable concern regarding our measurement is 
whether or not we have fitted to only the BAO fea- 
ture, or whether the result is offset due to noise spikes, 
which is obviously related to the significance of detec- 
tion. The conventional method of determining signifi- 
cance of BAO detection is to use a template with (i.e., 
BAO template) and without BAO information (i.e., no- 
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TABLE 2 

Best fit distance scale. 



z 


0.54 




Assumptions 


'"s,fld 


153.14 Mpc 




using Eisenstein & Hu figgS) 




[DA{z)lrs]iid 


8.584 


Results 


Best fit a 
Best fit DA(z)/rs 
Best fit Da(z) 


i.UDDl_(j g4gg 

1411 ± 65 Mpc 
using the prior, Ts = 153.7 ± 1.7 Mpc (WMAP7) 



Note. — Our best fit BAO scale and the derived distance scale. 




0.2 0.4 0.6 0.8 1 



z 



Fig. 9. — Interpreting o—l{%) = 6.61. Strictly speaking, the best 
fit value of a represents a constant ratio of the observed Da{z)/ts 
to the fiducial [D a{z) / r s\fid w6 assume in the template. The black 
solid and dashed lines (with a shade) show what the best fit and 
the 1 — (T error on a imply in this strict interpretation. However, 
although the redshift dependence of Da{z) we assume spans z ~ 
— 1, most of the galaxies are within z ~ 0.45 and 0.65 with a 
peak of the distribution near z ~ 0.5 — 0.55 (the shaded region in 
the bottom panel). Therefore, it is reasonable to consider that the 
best fit a[D a(z) / r s\iid represents DA{z)/rs near z ~ 0.5 — 0.55. 
We also show the best fit using different template cosmologies: 
using w = —0.8 and —1. The bottom panel shows the ratios of 
different Dyi(2)/rs: Dyi/rs/[D^(2:)/rs]fld- One sees that the three 
different templates have a very similar shape in Da{z){z) over z = 
0.45 — 0.65, once the absolute difference is absorbed into a. The 
three templates return virtually the same DA{z)/rs at z = 0.54. 
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8.5 9 9.5 10 
D,(z = 0.54)/r^ 

Fig. 10. — The best fit Da{z = 0.54)/r3 for various template 
cosmologies. We also show the result using a different parametriza- 
tion (labeled with 'Concordance ACDM (AlBO)' using AlBO) and 
the result with systematics correction (labeled with 'Concordance 
ACDM (SYS)'). The one in the square box and the dotted vertical 
line are our fiducial choice: DA{z)/rs = 9.212+q'43j at z = 0.54. 
The results are consistent within 1%. 

featu re in the power spectrum m ay be wiped out by noise 
fe.g.. iCabre fc Gaztafiagal[201lD . Our mock test, Figure 
[Bl shows that, in the presence of the sample variance 
that is the same as our data, we recover the true acous- 
tic scale within ^ 4.6% in 68.3% of the time. Obviously, 
we are not fitting to a BAO feature in the tails of the 
distribution. We therefore rephrase our detection level: 
our measurement is likely to recover the true BAO scale 
within 4.6% in 68.3% of cases, assuming that BAO exists. 



BAO template) and observe the difference between 
values of the two best fits. Unfortunately, such Ax^ is 
very mo del- dependent. Using AOBl, we derive Ax^ = 4, 
which can be conventionally interpreted as a 2a detec- 
tion of BAO. However, such detection level depends on 
the choice of paramctrization. 

We reconsider this issue of the d etection level. Various 
observations including W MAP7 (iKoiiiatsu et al.ll20lH) 
and galaxy surveys (e.g., iBlake et al.l l2011bf ). have al- 
ready shown that BAO feature exists. Given the signal- 
to-noise ratio level of our data, we are interested in how 
likely we have fitted to a BAO feature not a noise feature, 
rather than detecting the existence of BAO. Therefore, 
rather than fitting the power spectra with a no-BAO tem- 
plate, we shall fit many realizations of the power spectra 
with a BAO template and sec how often we derive the 
correct BAO scale. For a large sample variance, the BAO 



6.3. Effect of systematics 

A number of observational systematics can potentially 
contaminate the observed galaxy clustering: stellar con- 
tamination, seeing variations, sky brightness variations, 
extinction, and color offsets (Schlafly et al. 2010). How- 
ever, as long as the systematics do not introduce a pre- 
ferred scale similar to the BAO scale, i.e., if the system- 
atics only introduce a smooth component in the power 
spectrum up to a sample variance, our results would not 
depend on the contamination from systematics. Paper 
I moreover has shown that the effect of the survey sys- 
tematics are small. We therefore have not included the 
systematic corrections for our main result. In this sec- 
tion, however, we use power spectra that were corrected 
for the sys tematics using th e method introduced in Pa- 
per I (See iRoss et ani2011l for a similar method for the 
correlation function) and observe the effect of the sys- 
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tematics on the result. 

The method in Paper I assumes that the effect of sys- 
tematics is small and linear. In Fourier space, there- 
fore, we assume that the following equation holds for 
each wave band £: 



(15) 



where Sz^ii) is the observed galaxy density field at the 
zf^ redshift bin, 6gj is the true galaxy density field, and 
are the variation of systematics across the sky. We 
only include the dominant three systematics identified 
in Paper I: stellar contamination, seeing variations, and 
sky brightness variations. If wc assume that there is no 
intrinsic correlation between the systematics and the un- 
derlying large scale structure, i.e., < (5g,zi<5s,a(^) >= 0, 
we can solve for ez^^a using the measurements of galaxy 

power spectra (i.e., < Sz^Szj >) and the cross-power 
spectra between galaxies and the systematics (i.e., < 

^Zi^sa >), as presented in Paper I. The error on the band 
power is minimally propagated: the error is quadratically 
increased by the amount of the final correction, after tak- 
ing into account the number of wave modes. 

Figure [TT] and [T^ show the best fit results when we use 
the power spectra after systematics correction. We de- 
rive a -1 = 7.0l2ti-l\%. The reduced has slightly 
improved to be 1.09. The difference in between us- 
ing the BAO template and the no-BAO template has 
increased to 6.2 after systematics correction, from the 
previous 4.0 without systematics correction. The right 
panel of Figure [121 in comparison to Figure [8l shows 
that the systematics correction, while the effect is small, 
improves the fit on large scales (/ < 100). 

Figure [10] shows that the best fit value with the system- 
atics correction is consistent with the fit before the sys- 
tematics correction within 1% of a, demonstrating that 
the BAO fitting is fairly robust against the systematics 
effects. The overall improvement in the statistics after 
the systematics correction, such as on the reduced x^i 
motivates the usage of the method in Paper I for fu- 
ture surveys, which can be further improved with a more 
careful error propagation during the correction. 

7. DISCUSSIONS: COSMOLOGICAL 
IMPLICATIONS 

We combine our measurements of Da{z = 0.54) with 
recent spectroscopic BAO measurements. The spectro- 
scopic surveys report Dy(z) that contains both the in- 
formation along the line of sight, H{z), and the infor- 
mation on the transverse direction. Da- In Figure [T3l 
wc present the distancc-to-rcdshift relations of different 
BAO measurements in a 2-dimcnsional space of Da{z) 
and H[z). Our measurement of Da{z = 0.54) appears 
as the black horizontal line with the shaded region repre- 
senting the associated error. We also show the measure- 
ments of Dyi z = 0.2) /rg and Dv{z = 0.35) /r., from 
Percival et alT ((20TO ) for SDSS DR7 ([Abazaiian et al.l 
as red lines with r nagenta shades and Dv{z = 
0.6) /rs from iBlake eTall (|2011bD for the WiggleZ data 
over 0.2 < z < 1 as a green line with a light green shade. 
The black square points (along the dotted line) show the 
expected Da{z) and H at z = 0.2, 0.35, 0.54, and 0.6 



1.5 



0.5 



Owg - l-OfO (- 0.045, + 0.047) 
Xl, = 84-54 
. CMASSl 




u 




Fig. 11. — The best fit result after systematics correction using 
the combinations of CMASSl, CMASS2, CMASS3, and CMASS4 
and using AOBl. We derive a - 1 = 7.012t*-^J%, which is quite 
similar to the result before systematics correction. The black data 
points in the four panels with error bars show the measured C'l 
after systematics correction divided by a smooth fit at CMASSl, 
CMASS2, CMASS3, and CMASS4. The red lines represent the 
best fit C; and the red circles show the best fit band power Ci 
after the window function effect is considered. 
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Fig. 12. — Left panel: The surfaces along a after systcmatics correction for Figure [TT] when marginalized over other parameters (red 
line). Th e bl ue line shows the surface when the BAO is removed from the template. Right panel: a stacked Ci/Ci of the four panels 
of Figure [TT] after systematies correction. To better visualize the BAO feature we measured, we shift the wavenumbers of the four power 
spectra by £'A(-2mcdian)/i^A(-2 = 0.54), re-bin the combined band powers while inversely weighting by errors. The solid red line is the best 
fit for CMASS2, after the wavenumber is rescaled to z = 0.54. 



TABLE 3 

The derived cosmological parameters. 



^ Bl'ake 2011b 
H Percival 2010 



h 

w 



L>ODM 



oODM 



0.2912 (0.2917) ± 0.0270 
0.6884 (0.6892) ± 0.0392 
-1.0185 (-1.0337) ±0.1862 
0.7088 (0.7083) ± 0.2705 
Fixed at Qk = 0.0 



0.2939 (0.2952) ±0.0170 
0.6748 (0.6715) ±0.0175 

Fixed at to = —1.0 
0.7118 (0.7116) ±0.0172 
-0.0057 (-0.0067) ± 0.0058 



Note. — Marginalized fit and errors associated with the 
fit on selective param eters that are derived using COSMOMC 
l|Lewis fc Bridld 120021 ) for two different cosmologies. The value 
inside the parentheses show the best fit values. 



for our fiducial ACDM. Note that the measurements be- 
yond z = 0.35 have a tendency to imply the location of 
the BAO at a smaller scale than the concordance ACDM 
(i.e., a larger Da{z) than the fiducial cosmology), in- 
cluding our Da measurement (~ lAa away). Due to 
nonlinear structure formation and galaxy bias, we ex- 
pect about a ~ 0.5% of bias t owards a smaller value 
on t h e measured BAO scale (ICrocce fc Scoccimarrol 
200l IPadmanabhan fc Whit^ 120091 : ISeo et al.l 120101: 
Mehta et al.l l201l[) . which has not been accounted in 
these measurements. Such correction will slightly im- 
prove the consistency between the BAO measurements 
and the concordance ACDM, but it is overall a very small 
effect for the current level of errors. The circles along the 
dashed line and the crosses along the dot-dashed line in 
Figure [13] show the expected Da{z) and H{z) based on 
our best fit wCDM and oCD M cosmologies from COS- 
MOMC (jLewis fc Bridlell2002D that will be explained be- 
low. 

We use COSMOMC (ILewis fc Bridld[2002h to combine 
BAO measurements from the various gala xy surveys with 
the WMAP7 data (jKomatsu et al.ll201lD to derive con- 
straints on cosmological parameters. For BAO measure- 
ments, we use Dyiz = 0.2) /rs and Dv{z = 0.35)/rs 
from SDSS DR7 (|Percival et al.l[20lol) . Dyiz = 0.44)/rs, 



Q 




5 10 

cz/H(z)r 



Fig. 13. — Various BAO measurements in comparison to the 
concordance ACDM. The measurement of D^{z = 0.54) from 
this paper is shown with the black horizontal line. The gray 
shade represents I — a error. Red lines with magenta shades show 
Dy{z = 0.2) /rs and Dv(z = 0.35)/rs from Percival et al. (201(3) 
and the green line shows Dy(z = 0.6)/ra from Blake et al. (2011b|). 
The black squares along the diagonal dotted line show the expected 
combination of Da (z) and H based on the concordance ACDM at 
the redshifts of the data. One sees that the data beyond z = 0.35 
observed the BAO at a slightly smaller scale (i.e., a larger dis- 
tance) than the concordance ACDM. The circles along the dashed 
line and the crosses along the dot-dashed line show the expected 
Da{z) and H based on the best fit uiCDM and oCDM cosmologies 
in Table [3] 

Dviz = 0.60)/r,, and Dyiz = 0.73)/r, from WiggleZ^^, 
and DA{z){0.54)/rs from this work. The WMAP7 data 

For COSMOMC, we use the three - redshif t slice representation 
of the WiggleZ data from lBlake et al.l l l2011bl ), i.e., 0.2 < z < 0.6, 
0.4 < z < 0.8, and 0.6 < z < 1 .0, accounting for the covariance 
among them, while in Figure [131 we show the result for the whole 
redshift range (0.2 < z < 1.0). Note that the distance measure- 
ments from the WiggleZ data include non-BAO information. 
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provides the sound horizon scale and the distance to the 
last scattering surface and therefore, in combination of 
the BAO measurements from the galaxy surveys, we can 
break the degeneracies and constrain w and Vim (for 
luCDM) or 17a and VL,n (for oCDM). The cosmological 
parameters that the COSMOMC chain vary are ilbh^, 
^ch^, 6, T, Us, \nAs, and Asz, in addition to w (for 
wCDM) or riK (for oCDM); here, flch"^ is the dark mat- 
ter density, 9 is the approximate ratio of the sound hori- 
zon scale to the angular diameter distance to recombi- 
nation, r is the optical depth to reionization, As is the 
primordial superhorizon power in the curvature pertur- 
bation on 0.05 Mpc~^ scales, and Asz is the amplitude 
of the SZ power spectrum. 

The left panels of Figure [14] show marginalized 2- 
D likelihood contour surfaces that enclose 68.3% and 
95.5% of the likelihood (reddish shaded contours) on 
flm and w (top) and flm and h (bottom) assuming a 
flat wCDM, in comparison to the case without our mea- 
surement (dashed green lined contours for the spectro- 
scopic BAO measurements). The reddish contour lines 
in the top left show the constraint from our measure- 
ment alone using the current CMB prior on flmh^ (i-e, 
0.1326 ± 0.0063). This contour implies that, given the 
strong prior on Q^h^, adding our measurement of the 
distance scale at z = 0.54, which is larger than what 
is expected in the concordance ACDM, weighs toward a 
slightly larger 0^ and therefore a slightly smaller h with 
respect to the other data sets. We present the marginal- 
ized and the best fits of selective cosmological parameters 
in TableU = 0.2912 ± 0.0292, w = -1.0185 ± 0.186, 
h = 0.6884 ± 0.0392 for a flat wCDM. The right pan- 
els show the 2-D contour on il,„ and Ha for oCDM 
while holding w ~ —1. The best fit parameters are 
= 0.2939 ± 0.0170, = -0.0057 ± 0.0058, and 
h = 0.6748 ± 0.0175 in this case. Overall, an addition of 
our measurement slightly increases flm and decreases fix 
toward a more negative value. In terms of errors, includ- 
ing our data point provides only a slight improvement on 
ilm and h for oCDM. 

8. CONCLUSION 

We have measured the acoustic scale from the SDSS- 
III DR8 imaging catalog using 872,921 galaxies over 
~ 10,000deg^ between 0.45 < z < 0.65. Galaxies 
are binned into four different redshift slices where the 
width of each slice is 0.05, which is approximately the er- 
ror associated with photometric redshift determination. 
Angular power spectra are generated using an optimal 
quadratic estimator, as presented in Paper I. We use 
~ 110,000 SDSS III BOSS galaxies as a training sam- 
ple to derive the true redshift distribution of the galax- 
ies in the imaging catalog and therefore build reason- 
able template power spectra. We fit the templates to 
the measured angular power spectra and derive the best 
fit acoustic scale while marginalizing over sufficient free 
parameters to exclude any non-BAO signal. 

We derive DAiz)/rs = 9.212t|];^'j'^ at z = 0.54. Using 



the current WMAP7 constraint on the sound horizon at 
drag epoch, 153.2 ± 1.7 Mpc, we derive angular diameter 
distance -Da(z) = 1411 ± 65 Mpc at z = 0.54. Without a 
BAO feature in the template power spectrum, we cannot 
constrain a distance scale; the distance information we 
derive is therefore dominated by the BAO feature for our 
choice of parametrization. 

Our measurement of the distance scale is quite insen- 
sitive to the fiducial cosmology we assume for building 
the template. For a wide range of cosmologies we have 
investigated in this paper, the best fit varies less than 1% 
in the acoustic scale while the la error is ^ 4.7%. 

The angular distance scale we derive is 1.4cr higher 
than the concordance ACDM model. When combined 
with three other BAO measurements fro m SDSS DR7 
spect roscopic survey s at z = 0.2 and 0.3 5 (jPercival et al.l 
[20T0t) and WiggleZ ()Blake et al.ll2011b[) at z - 0.6, we 
find a tendency of cosmic distances measured using BAO 
to be larger than the concordance ACDM for z > 0.35. 
Adding our measurement with these BAO measurements 
in the presence of WMAP7 prior therefore shifts the best 
fit rt„i slightly larger than the concordance cosmology. 

In this paper, we have aimed at deriving a robust and 
conservative BAO information from the angular cluster- 
ing of galaxies. We find that an accurate determination 
of the true redshift distribution of galaxies is crucial for 
a good photometric BAO measurement. Although the 
details of the method would and should vary for the con- 
ditions of different surveys, we hope that the approach 
described in this paper serves as a valuable reference for 
the analyses of future photometric BAO surveys. 
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